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Abstract 

A  mathematical  model  predicting  the  performance  of  planar  solid  oxide  fuel  cell  (SOFC)  has  been  developed.  The  model  is  fuel  flexible,  which 
implies  not  only  pure  H2  but  also  any  reformate  composition  (H2,  H2O,  CO,  CO2  and  CH4)  can  be  used  as  a  fuel.  The  important  characteristic 
of  this  model  is  the  consideration  of  reaction  zone  layers  as  finite  volumes.  Reaction  zone  layers  are  thin  layers  in  the  vicinity  of  the  electrolyte 
where  electrochemical  reactions  takes  place  to  produce  electrons,  oxide  ions  and  water  vapor  (and/or  carbon  dioxide).  In  addition,  the  effect  of 
Knudsen  diffusion  is  accounted  in  the  porous  electrode  (backing)  and  reaction  zone  layers.  The  model  can  predict  the  performance  of  SOFC 
at  various  operating  and  design  conditions.  The  predicted  performance  of  SOFC  is  validated  with  the  measured  data  found  in  the  literature.  An 
excellent  agreement  is  obtained  between  the  predicted  performance  and  the  experimental  results.  Moreover,  the  effect  of  various  operating  and 
design  parameters  on  the  performance  of  SOFC  has  been  examined.  It  is  found  that  the  anode  concentration  overpotential  in  an  anode-supported 
SOFC  is  about  four  orders  of  magnitude  smaller  than  the  anode  ohmic  overpotential  even  at  higher  current  densities.  Further,  it  is  found  that  the 
single  largest  contributor  to  the  total  cell  potential  loss  is  ohmic  overpotential  and  hence  it  needs  to  be  minimized  to  enhance  the  cell  performance. 
©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Due  to  flexibility  in  fuel  choice,  solid  oxide  fuel  cells 
(SOFCs)  are  receiving  considerable  attention  for  both  small-  and 
large-scale  applications  [1],  Not  only  pure  Tb  but  also  any  refor¬ 
mate  composition  consists  of  multi-component  species  such  as 
H2,  H2O,  CO,  CO2  and  CH4  can  be  used  as  a  fuel.  Moreover, 
SOFCs  can  operate  directly  on  hydrocarbon  fuels  with  or  with¬ 
out  internal  reforming,  thereby  reducing  the  cost  of  an  external 
reformer  [2].  Hence,  fuel  flexibility  is  one  of  the  greatest  advan¬ 
tages  of  SOFCs  as  compared  to  other  types  of  fuel  cells. 

The  planar-type  design  of  SOFC  has  the  potential  to  offer 
higher  power  density  than  the  tubular  design  [3].  Due  to  its  com¬ 
pactness,  it  can  be  stacked  in  resemblance  to  polymer  electrolyte 
membrane  (PEM)  fuel  cells  to  satisfy  the  power  requirement 
of  an  application  [4],  Higher  power  density  of  planar  SOFC  is 
due  to  the  shorter  current  paths  resulting  in  low  ohmic  overpo¬ 
tential.  Moreover,  planar  SOFC  is  simple  to  fabricate  and  can 
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be  manufactured  into  various  configurations  [5].  However,  the 
problems  associated  with  the  development  and  commercializa¬ 
tion  of  planar  SOFCs  are  requirement  of  high  temperature  gas 
seals,  internal  stresses  in  cell  components  due  to  non-uniform 
temperature  distributions  and  high  manufacturing  cost. 

In  order  to  overcome  the  problems  associated  with  planar 
SOFCs,  much  of  the  recent  efforts  are  devoted  to  develop  new 
materials  and  configurations  to  improve  the  performance  at 
reduced  operating  temperatures  [2,6-9].  By  lowering  the  oper¬ 
ating  temperature  of  SOFC  to  around  700  °C,  many  of  the 
problems  associated  with  planar  SOFC  can  be  resolved.  For 
instance,  conventional  stainless  steel  can  be  used  for  intercon¬ 
nectors  instead  of  more  expensive  high  chrome  alloys  or  oxides 
and  hence  results  in  minimizing  the  material  and  manufactur¬ 
ing  costs  [10,11].  However,  the  ionic  conductivity  of  the  elec¬ 
trolyte  decreases  with  the  reduction  of  operating  temperature. 
The  reduction  in  ionic  conductivity  of  the  electrolyte  resulting 
in  higher  ohmic  overpotential  at  reduced  operating  temperatures 
can  be  minimized  by  either  using  electrode  (anode  or  cathode) 
supported  configuration  of  SOFC,  wherein  thin  electrolytes  of 
thicknesses  in  the  range  of  10-20  p.m  are  deposited  on  the  thick 
electrode  (anode  or  cathode)  or  using  composite  ceramic  elec- 
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Nomenclature 

Av 

reactive  surface  area  per  unit  volume  (m2  m-3) 

^*H2 

hydrogen  concentration  (molm-3) 

cH2,ref 

reference  hydrogen  concentration  (mol  m-3) 

OO2 

oxygen  concentration  (mol  m-3) 

c02,ref 

reference  oxygen  concentration  (mol  m-3) 

Dij 

binary  diffusion  coefficient  (m2  s- 1 ) 

OKn.i 

Knudsen  diffusion  coefficient  (m2  s-1) 

< 

effective  diffusion  coefficient  (m2  s-1) 

F 

Faraday’s  constant  (96,487  C  mol-1) 

J 

current  density  (A  m-2) 

Je 

electronic  current  density  (Am-2) 

Ji 

rtt2 

''O.ref 

ionic  current  density  (Am-2) 

reference  exchange  current  density  for  H2 
oxidation  (A  m-2) 

t02 

J0,ref 

reference  exchange  current  density  for  O2 
reduction  (Am-2) 

^bs 

backward  reaction  rate  constant  for  shift  reaction 
(molm-3  Pa-2  s-1) 

kfs 

forward  reaction  rate  constant  for  shift  reaction 
(molm-3  Pa-2  s-1) 

Kps 

equilibrium  constant  for  shift  reaction 

Mi 

molecular  weight  of  species  i  (kg  mol-1) 

n 

moles  of  electrons  transferred  per  mole  reactant 

nt  1 

number  fraction  of  electron  conducting  particles 
in  the  reaction  zone  layers 

«io 

number  fraction  of  ion  conducting  particles  in  the 
reaction  zone  layers 

nt 

total  number  of  particles  in  the  reaction  zone 
layers  per  unit  volume  (m-3) 

Ni 

diffusive  flux  of  species  i  (molm-2  s-1) 

P 

pressure  (Pa) 

Pel 

probability  of  percolation  of  electron  conducting 
particles  in  reaction  zone  layers 

Pio 

probability  of  percolation  of  ion  conducting 
particles  in  reaction  zone  layers 

A  1 

radius  of  electron  conducting 
particles  in  the  reaction  zone  layers  (m) 

^io 

radius  of  ion  conducting  particles  in  the  reaction 
zone  layers  (m) 

rs 

volumetric  shift  reaction  rate  (kgm-3  s-1) 

R 

universal  gas  constant  (8.3143  J  mol-1  K-1) 

volumetric  current  density  produced  due  to  H2 
oxidation  (Am-3) 

1ZC 

volumetric  current  density  produced  due  to  O2 
reduction  (A  m-3) 

Ss,i 

rate  of  production  or  consumption  of  species  i 
(kg  m-3  s- 1 ) 

ta 

anode  electrode  (backing)  thickness  (m) 

tarz 

anode  reaction  zone  thickness  (m) 

cathode  electrode  (backing)  thickness  (m) 

tcrz 

cathode  reaction  zone  thickness  (m) 

electrolyte  thickness  (m) 

T 

temperature  (K) 

Top  operating  temperature  (K) 

Vi  diffusion  velocity  of  species  i  (ms-1) 

x  rectangular  coordinate  (m) 

x,  mole  fraction  of  species  i 

Z  average  coordination  number 

Ze i  coordination  number  of  electron  conducting 

particles  in  the  reaction  zone  layers 
Zi0  coordination  number  of  ion  conducting  particles 
in  the  reaction  zone  layers 

Greek  symbols 

a  charge  transfer  coefficient 

yu2  reaction  order  for  Ht  oxidation 

yo2  reaction  order  for  O2  reduction 

e  porosity 

rj  overpotential  (V) 

k  ionic  conductivity  (S  m- 1 ) 

Pi  partial  density  of  species  i  (kg  m-3) 

a  electronic  conductivity  (S  m- 1 ) 

r  tortuosity 

<j)c  electronic  potential  (V) 

(pi  ionic  potential  (V) 

0  volume  fraction  of  electron  conducting  particles 
in  the  reaction  zone  layers 

Subscripts 

a  anode  activation 

bl  electrode  (backing)  layer 

c  cathode  activation 

rl  reaction  zone  layer 

Superscripts 
eff  effective 

e-  electron  conducting  particles 

O2-  ion  conducting  particles 


trolyte  in  an  electrolyte  supported  configuration  having  high 
ionic  conductivity  at  reduced  temperatures  [2]. 

Mathematical  modeling  is  an  essential  aspect  of  SOFC  tech¬ 
nology  development  process.  A  numerical  model  facilitates 
research  and  development  by  minimizing  the  need  of  repeti¬ 
tive  and  costly  experimentation.  Numerical  simulation  of  SOFCs 
provides  a  thorough  understanding  of  how  cell  performance  is 
affected  by  various  operating  and  design  parameters  such  as 
temperature,  pressure,  fuel  composition,  reactive  surface  area, 
pore  size,  porosity,  tortuosity,  thickness  of  various  components, 
and  ionic  and  electronic  conductivities  of  ion-conducting  and 
electronic-conducting  particles  of  the  electrode,  respectively, 
and  thus  helps  in  optimizing  cell  and  stack  design. 

Numerous  models  of  SOFCs  exists  in  the  literature 
[3-5,7,12-14],  varying  in  the  number  of  assumptions  employed. 
The  common  assumption  in  the  existing  SOFC  models  is  the 
consideration  of  reaction  zone  layers  as  mathematical  surfaces, 
treating  them  as  boundary  conditions.  However,  for  compos¬ 
ite  electrodes  such  as  those  in  SOFCs,  the  reaction  zone  layers 
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are  spreaded  out  into  the  electrode  some  distance  (10-50  p,m) 
from  the  electrolyte/electrode  interfaces  [15-17].  Reaction  zone 
layers  are  relatively  thin  layers  where  fuel  and  oxidant  are  elec- 
trochemically  converted  into  electrical  work,  heat  and  water 
vapor  (and/or  carbon  dioxide).  Although  Kim  et  al.  [7]  and  Zhu 
and  Kee  [  14]  discussed  about  reaction  zone  layers  but  no  one 
attempted  to  model  these  regions  as  finite  volumes,  rather  treated 
them  as  boundary  conditions. 

Of  the  models  existed  in  the  literature,  only  Kim  et  al.  [7] 
and  Zhu  and  Kee  [14]  validated  their  models  with  the  mea¬ 
sured  data  available  in  the  literature.  The  model  of  Kim  et  al. 
[7]  is  an  empirical  model  with  five  fitting  parameters  and  caters 
only  binary  mixtures  in  the  electrodes.  In  addition,  none  of  the 
existing  complete  cell  models  on  SOFC  incorporate  micro  char¬ 
acteristics  of  electrodes.  Incorporation  of  micro  parameters  into 
the  mathematical  cell  model  not  only  helps  in  improving  the 
accuracy  of  the  model  but  also  helps  in  quantitative  design  and 
optimization  process  [18]. 

The  objective  of  the  present  study  is  to  develop  a  mathemat¬ 
ical  model  of  SOFC  satisfying  the  following  requirements:  (i) 
flexibility  in  fuel  choice,  not  only  pure  H2  but  also  any  reformate 
composition  (FF,  H2O,  CO,  CO2  and  CH4)  can  be  used  as  a  fuel, 
(ii)  finite  reaction  zone  layers,  and  (iii)  incorporation  of  the  micro 
characteristics  of  the  electrodes.  Moreover,  the  other  objectives 
are  to  validate  the  developed  model  with  the  measured  perfor¬ 
mance  data  available  in  the  literature  and  to  examine  the  effect 
of  various  operating  and  design  parameters  on  the  performance 
of  SOFC. 

2.  Model  formulation 

Fig.  1  illustrates  different  layers  of  SOFC.  A  typical  SOFC 
consists  of  three  major  layers  such  as  anode  electrode  layer,  elec¬ 
trolyte  layer  and  cathode  electrode  layer.  Two  additional  layers, 
which  are  shown  in  Fig.  1  between  the  anode  electrode  layer 
and  electrolyte  layer  and  cathode  electrode  layer  and  electrolyte 
layer  are  called  anode  reaction  zone  layer  and  cathode  reaction 
zone  layer,  respectively.  The  reaction  zones  on  either  side  of 
the  electrolyte  are  considered  to  be  distinct  layers  because  of 
the  fact  that  the  electrochemical  reaction  not  only  occurs  on  the 
interface  between  the  anode  and  electrolyte,  and  the  cathode 
and  electrolyte  but  also  extends  to  a  depth  of  10-50  pm  inside 
electrode  layers  [15-17]. 

The  cell  is  assumed  to  operate  under  steady  state  condition 
and  the  parameters  vary  in  the  .v-direction  only,  as  shown  in 
Fig.  1.  The  temperature  and  total  pressure  are  assumed  to  be 
uniform  in  the  electrode  and  reaction  zone  layers.  The  convec¬ 
tive  flux  is  assumed  to  be  negligible  in  the  porous  electrode 
and  reaction  zone  layers  when  compared  to  the  diffusive  flux  of 
gaseous  species,  which  implies  that  the  primary  means  of  species 
transport  in  the  porous  electrode  and  reaction  zone  layers  is  dif¬ 
fusion  [7].  The  reactant  gas  mixtures  are  approximated  as  ideal 
gases  with  negligible  viscous,  Soret,  Dofour  and  gravity  effects. 
Since  the  reaction  zone  layers  are  considered  as  separate  regions, 
there  are  no  electrochemical  reactions  (either  oxidation  or  reduc¬ 
tion)  in  the  electrode  layers.  The  reaction  zone  layers  consists  of 
mixture  of  electron-conducting  particles,  electrolyte-conducting 
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Fig.  1.  Illustration  of  different  layers  of  solid  oxide  fuel  cell. 

particles  and  void  space  occupied  by  gaseous  species.  Moreover, 
in  the  present  study,  it  is  assumed  that  there  is  no  electrochemical 
oxidation  of  CO  in  the  anode  reaction  zone  layer.  The  electrolyte 
layer  is  assumed  to  be  a  dense  solid  with  no  interconnected 
porosity.  With  these  assumptions  the  cell  model  is  formulated 
and  is  described  in  the  following  sections. 

2.1.  Electrode  (backing)  layers 

The  processes  that  need  to  be  modeled  in  the  electrode  lay¬ 
ers  are  transport  of  multi-component  species  to  and  from  the 
reaction  sites  in  the  reaction  zone  layers  and  electron  migration 
through  the  solid  portion  of  the  porous  structure.  In  Fig.  1,  the 
electrode  layers  correspond  to  regions  a-b  and  e-f. 

The  mathematical  model  for  the  processes  occurring  in  the 
electrode  layers  are  formulated  by  applying  the  conservation  of 
species  and  conservation  of  electronic  charge.  The  governing 
equations  are  described  as  follows: 


species  :  V  ■  (p,V/)  =  Ss.i 

(1) 

electronic  charge  :  V  •  /e  =  0 

(2) 

where  p,-  is  the  partial  density  of  species  i,  V, ■  is  the  species 
diffusion  velocity,  is  the  species  source  term,  and  Je  is  the 
electronic  current  density. 

The  diffusion  velocity  or  in  other  words,  the  diffusion  flux  in 
the  SOFC  electrode  layers  can  be  determined  using  the  modified 
Stefan-Maxwell  equations  incorporating  Knudsen  diffusion  for 
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multi-component  systems  involving  n  species  [19]: 

n  | 

-c  Vx;  =  (xjN<  ~  xiNj )  (3) 


'  D 

j=U  >J 


where  c  is  the  concentration  of  the  mixture,  x,  the  mole  fraction  of 
species  i,  N,  the  diffusion  flux  of  species  i,  and  Df/  the  effective 
diffusion  coefficient,  defined  as: 


Dcfi  _  ®  f  DijDKn,i  \ 

U  T  VAy  +  OKn,// 


(4) 


where  Du  is  ordinary  or  binary  diffusion  coefficient  and  Z>Kn,/  is 
Knudsen  diffusion  coefficient  of  species  i,  e  and  r  are  porosity 
and  tortuosity  of  the  electrode  layers,  respectively. 

The  species  source  term  on  the  right-hand  side  of  Eq.  (1) 
represents  the  rate  of  production  or  consumption  of  species  due 
to  chemical  reactions.  Since  there  is  no  chemical  reaction  in 
the  cathode  electrode  layer,  the  species  source  term  in  Eq.  (1) 
vanishes,  whereas  in  the  anode  electrode  layer,  it  represents  the 
volumetric  production  or  consumption  of  species  due  to  water- 
gas  shift  reaction  given  as: 


CO  +  H20  ^  H2  +  C02 


(5) 


The  volumetric  reaction  rate  for  the  shift  reaction  can  be 
written  as: 


rs  =  A:fsccoCH20  -  ^bscco2cH2 


(6) 


where  kfs  and  k\,s  are  forward  and  backward  reaction  rate  con¬ 
stants  for  the  water-gas  shift  reaction. 

In  terms  of  mole  fraction,  the  volumetric  reaction  rate  for 
shift  reaction  becomes: 

rs  —  (^t)  [£fs-*C0*H20  _  &bs*C02*H2]  (7) 

where  p  is  the  total  pressure,  R  is  the  universal  gas  constant  and 
T  is  the  temperature. 

The  mass  rates  of  production  or  consumption  of  various 
species  in  the  anode  electrode  layer  can  be  formulated  as  fol¬ 
lows: 


Ss,h2  =  A%Ts  (8) 

‘S's.HjO  =  (9) 

Ss.co  =  -Mco^s  (10) 

Ss,co2  =  MCQ2rs  (11) 


where  Mh2,  Mh2o>  Mqo  ,  andMco2  are  the  molecular  weights 
of  H2,  H20,  CO,  and  C02,  respectively. 

Using  Ohm’s  law,  the  electronic  current  density  given  in  Eq. 
(2)  can  be  written  in  terms  of  electronic  potential  in  the  electrode 
layers  as: 

electronic  charge  :  V  ■  (oyffV0e)  =  0  (12) 

where  aff  is  the  effective  electronic  conductivity  of  the  electrode 
layers,  defined  as: 


eff  . 


1  —  E 

r 


a 


(13) 


where  a  is  the  electronic  conductivity  of  pure  electrode  layer 
materials. 


2.2.  Reaction  zone  layers 

As  mentioned  in  the  previous  sections,  reaction  zone  layers 
are  relatively  thin  layers  where  fuel  and  oxidant  are  electrochem- 
ically  converted  into  electrical  work,  heat  and  water  vapor.  Often, 
these  layers  are  idealized  as  mathematical  surfaces,  where  the 
electrochemical  reactions  (oxidation  or  reduction)  are  assumed 
to  occur  at  the  localized  triple  phase  boundary  (TPB)  between 
the  electron-conducting  particle,  the  ion-conducting  particle  and 
the  reactant  gas  atmosphere.  However,  consideration  of  reac¬ 
tion  zone  layers  as  finite  volumes  would  increase  the  number 
of  reactive  sites  and  leads  to  better  agreement  between  the 
model  prediction  and  the  experimental  observation  [20].  There¬ 
fore,  in  this  study  the  reaction  zone  layers  are  considered  as 
finite  volumes  rather  than  mathematical  surfaces,  comprising  of 
electron-conducting  particles,  ion-conducting  particles  and  void 
spaces  occupied  by  the  gases.  The  processes  which  need  to  be 
modeled  in  these  layers  are  transport  of  gaseous  species  with 
chemical  and  electrochemical  reactions  and  electron  and  ion 
migration  in  the  electron-  and  ion-conducting  particles,  respec¬ 
tively. 

The  conservation  of  species  equation  governing  the  trans¬ 
port  of  multi-component  species  in  the  reaction  zone  layers 
is  similar  to  the  electrode  layers  given  in  Eq.  (1)  except  the 
source  term,  which  also  accounts  for  species  consumption  or 
production  due  to  electrochemical  reactions.  Again,  the  modified 
Stefan-Maxwell  equation  given  in  Eq.  (3)  is  used  to  model  dif¬ 
fusion  flux  in  the  reaction  zone  layers.  The  equations  governing 
the  migration  of  electrons  and  ions  in  the  electron-  and  ion¬ 
conducting  particles  of  the  reaction  zone  layers  are  described  as 
follows: 


•  anode  reaction  zone  layer: 

electronic  charge  :  V  •  Je  =  TZa  (14) 

ionic  charge  :  V  •  7;  =  —7 Zd  (15) 

•  cathode  reaction  zone  layer: 

electronic  charge  :  V  ■  Je  =  —1ZC  (16) 

ionic  charge  :  V  •  j;  =  1ZC  (17) 


where  lZd  and  1ZC  are  volumetric  current  densities  produced  in 
the  anode  and  cathode  reaction  zone  layers  due  to  H2  oxidation 
and  02  reduction  reactions,  respectively,  which  are  given  by  the 
general  Butler-Volmer  equations  as: 


lZd  —  Av/^ef 


CH2\ra2 

CH2,ref/ 


X 


an  Fijd 
RT 


(1  -  a)nFij.d  ^  | 


(18) 
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TZC  —  AyjO?ef 


exp 


^  ( )■ 
c02.ref 
cm  Frjc 
RT 


YO-, 


-exp  - 


Cl  -  a)nFrjc 
~RT 


Pio  = 


(19)  where 


where  Av  is  the  actual  reactive  surface  area  per  unit  volume, 
and  /q  ,j?ef  are  reference  exchange  current  densities  for  H2  oxida¬ 
tion  and  O2  reduction  at  the  reference  concentrations  CH-,,ref  and 
co2,ref,  respectively,  a  the  charge  transfer  coefficient  (or  sym¬ 
metry  factor),  whose  value  lies  between  0  and  1,  n  the  number 
of  electrons  participating  in  the  electrochemical  reaction  and  ija 
and  tjc  are  anode  and  cathode  activation  overpotentials,  respec¬ 
tively,  defined  as: 


0a  —  $1  0e 
0c  —  0e 


(20) 

(21) 


The  expression  used  to  model  reactive  surface  area  per  unit 
volume  (Av)  was  developed  by  Costamagna  et  al.  [21]  based  on 
the  particle  coordination  number  in  binary  random  packing  of 
spheres,  given  as: 


X  •  2  /I  T  Zel  Z  io 

Av  =  jt  sin  9relntnelnio — —  peiPio 


(22) 


where  6  is  the  contact  angle  between  electron  and  ion  conducting 
particles  in  the  reaction  zone  layer,  rei  is  the  radius  of  the  electron 
conducting  particles,  nt  is  the  total  number  of  particles  per  unit 
volume,  ne\  and  nm  are  the  number  fraction  of  electron  and  ion 
conducting  particles  in  the  reaction  zone  layers,  respectively, 
Ze  1  and  •Zio  are  the  coordination  numbers  of  electron  and  ion 
conducting  particles  in  the  reaction  zone  layers,  respectively, 
Z  is  the  total  average  number  of  contacts  of  each  particle  and 
Pei  and  p\„  are  the  probabilities  of  electron  and  ion  conducting 
particles  in  the  reaction  zone  layers,  respectively. 

The  parameters  required  to  obtain  reactive  surface  area  per 
unit  volume  (Av)  can  be  calculated  from  [18,21]: 


1  —  s 


n  t  = 


(4/3)7Tr|j [«el  +  (1  -  «el)(no/rel)3] 
where  e  is  the  porosity  of  the  reaction  zone  layers,  and: 
<P 

«e  1 


[0+((l-0)/(rio/rei)3)] 


(23) 


(24) 


where  <t>  is  the  volume  fraction  of  the  electron  conducting  par¬ 
ticles  in  the  reaction  zone  layers. 


Zel  —  3  + 


Z-  3 


Z,o  =  3  ■ 


[«e  1  +  (1  -  Hei)(rio/ r el)“] 

(Z  -  3)(ri0/rei)2 
[nei  +  (1  -  nei)(rio/rei)2] 


(25) 


(26) 


where  Z  is  the  total  average  coordination  number,  equal  to  6 
[18]. 


1-  2- 


-Zel— el 


2.5 


0.4 


(27) 


2.5 


-1  0.4 


■Zel— el  — 


«el  Z2 


Zin— in  — 


HinZi 


10  io 


(28) 


(29) 


(30) 


The  species  source  terms  in  the  anode  and  cathode  reaction 
zone  layers  due  to  electrochemical  reactions  are  related  to  volu¬ 
metric  current  density  through  Faraday’s  law  of  electrochemical 
reaction  and  are  described  as  follows: 


•  anode  reaction  zone  layer: 
ViMiTla 


$s,i  — 


nF 


(31) 


where  M;  is  the  molecular  weight  of  the  species  i,  7Zd  is  the  vol¬ 
umetric  current  density  produced  due  to  H2  oxidation  reaction 
given  in  Eq.  (18),  n  is  the  number  of  electrons  participating  in 
the  electrochemical  reaction,  F  is  the  Faraday’s  constant,  and 
Vj  is  the  stoichiometric  coefficient  of  the  species  involved  in 
the  oxidation  reaction  expressed  in  the  following  form  [22]: 


H2  +  O2-  -  H20  -*  2e“ 

•  cathode  reaction  zone  layer: 
ViMiTZc 


Ss.i  — 


nF 


(32) 


(33) 


where  1ZC  is  the  volumetric  current  density  produced  due  to 
O2  reduction  reaction  given  in  Eq.  (19),  expressed  in  the  fol¬ 
lowing  form  [22] : 


^02-0 


2- 


2e“ 


(34) 


Converting  the  electronic  and  ionic  current  densities  given 
in  Eqs.  (14)— (17)  into  electronic  and  ionic  potentials  through 
Ohm’s  law,  Eqs.  (14)— (17)  becomes: 


•  anode  reaction  zone  layer: 


electronic  charge  :  V  •  (cr^Vt^e)  =  —7 Z- 


ionic  charge  :  V  ■  (<7yfV</)j)  =  7 Z: 


cathode  reaction  zone  layer: 


electronic  charge  :  V  •  (cr^Vt/ie)  =  7 Z, 


ionic  charge  :  V  •  (erjr^  V^;)  =  —7 Z, 


(35) 

(36) 

(37) 

(38) 


where  the  effective  electronic  and  ionic  conductivities  in  the 
reaction  zone  layers  are  defined  as  [23]: 


erf  =  0 


1 


Pel  = 


a 


(39) 
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=  (1  -  <p)  (^)  K  (4°) 

where  <P  is  the  volume  fraction  of  the  electron  conducting  par¬ 
ticles  in  the  reaction  zone  layer,  a  and  k  are  the  conductivities 
of  pure  electron  and  ion-conducting  materials,  respectively. 

2.3.  Electrolyte  layer 

Electrolyte  layer  in  SOFC  is  fully  dense  with  no  intercon¬ 
nected  porosity.  The  major  function  of  the  electrolyte  is  to 
conduct  oxide  ions  produced  in  the  cathode  reaction  zone  layer 
to  the  anode  reaction  zone  layer,  thus  completing  the  electrical 
circuit.  The  conservation  equation  governing  the  migration  of 
oxide  ion  in  the  electrolyte  layer  is  given  as: 

V  •  Ji  =  0  (41) 

where  J\  is  the  ionic  current  density,  which  is  equal  to  the  total 
current  density  J. 

Using  Ohm’s  law,  the  ionic  current  density  can  be  expressed 
in  terms  of  ionic  potential  and  the  above  conservation  equation 
becomes: 

V  •  OrV</>i)  =  0  (42) 

where  k  is  the  ionic  conductivity  of  the  electrolyte. 

The  governing  equations  used  in  the  model  formulation  of 
different  components  of  SOFC  are  summarized  in  Table  1 . 


2.4.  Boundary  conditions 


The  location  at  which  boundary  conditions  are  needed  to 
complete  the  mathematical  formulation  are  illustrated  in  Fig.  1 
as  “a”,  “b”,  “c”,  “d”,  “e”  and  “f”.  Locations  “a”  and  “f”  are  the 
interfaces  between  the  fuel  channel  and  anode  electrode  (back¬ 
ing)  layer  and  air  channel  and  cathode  electrode  (backing)  layer, 
respectively;  the  boundary  conditions  at  these  location  are  spec¬ 
ified  boundary  conditions,  where  the  composition  of  gaseous 
species  and  electronic  potential  are  specified.  The  boundary  con¬ 
ditions  at  the  interfaces  between  the  anode  electrode  and  anode 
reaction  zone  layers  (location  “b”)  and  the  cathode  electrode 


and  cathode  reaction  zone  layers  (location  “e”)  are  continu¬ 
ous  flux  boundary  conditions  for  gaseous  species  and  electronic 
potential,  where  the  diffusion  flux  and  electronic  current  density 
are  continuous,  whereas  insulated  boundary  condition  for  ionic 
potential,  which  implies  ionic  current  density  is  zero.  At  loca¬ 
tions  “c”  and  “d”,  the  interfaces  between  the  anode  reaction  zone 
layer  and  electrolyte  layer,  and  the  cathode  reaction  zone  layer 
and  electrolyte  layer,  the  diffusion  flux  and  electronic  current 
density  are  zero,  whereas  the  ionic  current  density  is  continu¬ 
ous  and  is  equal  to  the  total  current  density.  Mathematically, 
the  boundary  conditions  at  different  locations  in  Fig.  1  can  be 
expressed  as: 

at  x  =  “a”  and  “f”  :  x,  =  specified,  </>e  =  specified; 

at  x  —  b  and  e  :  N ) 1 =  A)j rz '  J& I h i  —  7e|ri,  Ji  —  0; 

at  x  —  “c”  and  “d”  :  Nt  =  0,  Je  =  0,  /;  =  J. 

3.  Model  validation 

Since  the  mathematical  model  presented  above  is  based  on  the 
conservation  principle,  the  numerical  method  used  to  discretize 
the  governing  equations  is  a  well-known  method,  often  referred 
to  as  “finite-volume  method”  or  “control-volume  method”.  The 
discretized  governing  equations  form  a  system  of  algebraic  equa¬ 
tions  whose  solution  can  be  obtained  using  a  direct  method  or  an 
iterative  method.  For  the  present  problem,  an  iterative  method 
is  used  because  of  the  fact  that  the  coefficients  of  the  discretized 
algebraic  equations  are  not  constants  but  are  functions  of  the 
variables  evaluated  at  grid  points. 

Using  the  parameters  listed  in  Table  2,  the  predicted  cell  per¬ 
formance  is  compared  with  the  experimental  data  found  in  the 
literature  [24],  and  is  shown  in  Fig.  2.  The  performance  of  the 
cell  is  predicted  when  the  cell  is  supplied  with  95%  Hi  and  5% 
HiO  as  fuel,  operating  at  a  temperature  and  pressure  of  1073  K 
and  1  atm,  respectively.  Oxygen  composition  in  the  ambient  air 
is  used  as  oxidant.  It  can  be  seen  that  the  present  model  pre¬ 
dictions  agree  well  with  the  experimental  results.  It  is  worth 
mentioning  that  all  the  parameters  used  in  the  model  validation 
are  obtained  from  Rogers  et  al.  [24]  except  the  parameters  given 
with  the  references.  The  value  of  tortuosity  is  varied  to  obtain  the 


Table  1 

Governing  equations 

Cell  layers 

Conservation  of  species 

Conservation  of  electronic  charge 

Conservation  of  ionic  charge 

Anode  electrode  layer 

V  ■  (NiMi)  =  SsJ 

n 

,  -Vx;  =  Y ~]jai(XjNi  -  XjNj) 

11 

v  •  (crbf  V0e)  =  0 

- 

Cathode  electrode  layer 

V  ■  (NiMi)  =  0, 

n 

-Vxi  =  y^ixjNi  -  XiNj) 
i=U  u 

V  •  (abf  V0e)  =  0 

- 

Anode  reaction  zone  layer 

V  ■  (NMi)  =  Ss,i , 

n 

•  ~vx> = ~  *'A) 

j= i, ij 

V  ■  (of' V0e)  =  -TZ, 

v  •  (ojf  V<h)  =  Ra 

Cathode  reaction  zone  layer 

V  ■  (NMi)  =  SSii 

n 

,  -Vx;  =  V'-i (XjNj  -  XiNj ) 

.  y 

V  ■  (crfVfc)  =  Tlc 

V  ■  (afVfr)  =  -nc 

Electrolyte  layer 

- 

J=U 

- 

V  ■  (*-V0i)  =  0 
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Table  2 

Parameters  used  for  model  validation 


Operating  temperature,  Top  (K) 

1073.0 

Total  pressure,  p  (atm) 

1.0 

Fuel  composition,  xu2 ;  xh2o 

0.95;  0.05 

Air  composition,  xo2  5  *n2 

0.21;  0.79 

Anode  conductivity,  o  (S  m-1) 

71428.57 

Cathode  conductivity,  o  (S  m-1) 

5376.34 

Electrolyte  conductivity,  k  (S  m-1) 

0.64 

Anode  electrode  layer  thickness,  ta  (|xm) 

1000.0 

Cathode  electrode  layer  thickness,  tc  (|xm) 

50.0 

Anode  reaction  zone  layer  thickness,  tarz  (|xm) 

20.0 

Cathode  reaction  zone  layer  thickness,  tCIZ  (|xm) 

20.0 

Electrolyte  thickness,  te  (|xm) 

10.0 

Porosity  of  anode  and  cathode,  e 

0.375 

Tortuosity  of  anode  and  cathode,  r 

2.75  [20] 

Pore  diameter  of  anode  and  cathode,  dp  (|xm) 

1.5  [7] 

Contact  angle  between  e-  and  O2- conducting 

15  [21] 

particles,  6  (°) 

Radius  of  e-  conducting  particles,  rei  (|xm) 

0.1  [18,29] 

Radius  of  O2-  conducting  particles,  ri0  (|xm) 

0.1  [18,29] 

Volume  fraction  of  e-  conducting  particles,  <P 

0.5  [21,23] 

Reference  H2  concentration,  <?h2  (mol  m-3) 

10.78 

Reference  O2  concentration,  co2  (mol  m-3) 

2.38 

Reference  exchange  current  density  for  H2 

1320  [7] 

oxidation,  7?^  (Am-2) 

Reference  exchange  current  density  for  O2 

400  [21] 

reduction,  Jq^  (Am-2) 

Reaction  order  for  H2  oxidation,  yr2 

0.5 

Reaction  order  for  O2  reduction,  yo2 

0.5 

Source :  Rogers  et  al.  [24]. 


best  agreement  between  the  present  model  predictions  and  the 
experimental  results  shown  in  Fig.  2,  since  the  value  reported 
by  Rogers  et  al.  [24]  is  unity.  The  typical  tortuosity  value  for 
SOFC  electrodes  is  in  the  range  of  2-6  [20,25].  Hence,  the  tor¬ 
tuosity  value  of  2.75  used  in  the  present  model  predictions  is 
in  the  typical  range  for  SOFC  electrodes  and  provides  the  best 
agreement  with  the  experimental  results.  The  A1’ -value,  an  indi¬ 
cator  between  0  and  1 ,  reveals  how  closely  the  predicted  values 


Fig.  2.  Comparison  between  the  model  predictions  and  experimental  results  of 
Rogers  et  al.  [24]. 


correspond  to  the  experimental  data,  is  obtained  for  the  present 
model  as  0.988. 

4.  Results  and  discussion 

The  developed  SOFC  model  can  be  used  to  investigate  the 
effect  of  operating  and  design  parameters  on  the  cell  perfor¬ 
mance.  The  fuel  composition  used  in  the  simulation  is  listed 
in  Table  3.  Oxygen  composition  in  the  ambient  air  is  used 
as  oxidant.  The  model  base  case  parameters  derived  from  the 
data  available  in  the  literature  are  given  in  Table  4.  The  most 
significant  difference  between  the  base  case  and  model  valida¬ 
tion  parameters  given  in  Table  2  is  the  thickness  of  the  anode 
electrode  layer.  Typical  thickness  of  the  anode  electrode  layer 
reported  in  the  literature  for  an  anode-supported  SOFC  is  2  mm 
[10,16,26].  Since  SOFCs  are  fuel  flexible,  methane  reformate 
fuel  composition  given  in  Table  3  is  used  in  the  present  simula¬ 
tion.  Further,  water-gas  shift  reaction  is  considered  in  the  anode 
electrode  and  reaction  zone  layers. 

The  cell  performance  resulting  from  the  base  case  parame¬ 
ters  is  shown  in  Fig.  3.  It  can  be  seen  that  the  most  significant 
contribution  to  the  cell  potential  loss  is  from  the  anode  side  of 
the  cell.  The  anode  overpotential  includes  the  anode  activation 
overpotential  due  to  the  resistance  to  the  charge  transfer  reaction, 
anode  ohmic  overpotential  due  to  the  resistance  to  the  flow  of 
electrons  and  ions  in  the  anode  reaction  zone  layer  and  resistance 
to  the  flow  of  electrons  in  the  anode  electrode  (backing)  layer 
and  anode  concentration  overpotential  due  to  the  resistance  to 


Table  3 

Mole  fraction  of  methane  reformate  used  in  the  simulation 


Mole  fraction 

Value 

-^Hl2 

0.410 

-»'h2o 

0.435 

*co 

0.090 

*co2 

0.065 

Source :  Rajesh  et  al.  [31]. 
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Table  4 

Base  case  parameters  used  in  the  simulation 


Operating  temperature,  Top  (K) 
Total  pressure,  p  (atm) 

Anode  electrode  layer  thickness, 
k  (p-m) 

Cathode  electrode  layer 
thickness,  tc  ([xm) 

Anode  reaction  zone  layer 
thickness,  ?arz  (|xm) 

Cathode  reaction  zone  layer 
thickness,  tc rz  (jxm) 

Electrolyte  thickness,  tt  (|xm) 
Anode  conductivity,  o  (S  m-1) 
Cathode  conductivity,  a  (S  m-1) 
Electrolyte  conductivity,  k 
(S  m_1) 

Porosity  of  anode  and  cathode,  s 
Tortuosity  of  anode  and  cathode, 
r 

Pore  diameter  of  anode  and 
cathode,  dp  (|xm) 

Contact  angle  between  e_  and 
O2-  conducting  particles,  6  (°) 
Radius  of  e_  conducting 
particles,  rt\  (|xm) 

Radius  of  O2-  conducting 
particles,  ri0  (|xm) 

Volume  fraction  of  e_  conducting 
particles,  <P 

Reference  H2  concentration,  cu2 
(molm~3) 

Reference  O2  concentration,  co2 
(mol  m~3) 

Reference  exchange  current 
density  for  H2  oxidation,  J^ef 
(Am"2) 

Reference  exchange  current 
density  for  O2  reduction,  /Pjef 
(Ami 

Reaction  order  for  Ht  oxidation, 

m2 

Reaction  order  for  O2  reduction, 

Yo2 

Forward  reaction  rate  constant 
for  the  shift  reaction,  kfs 
(mol  m~3  Pa-2  s_1) 

Backward  reaction  rate  constant 
for  the  shift  reaction,  k\,s 


1073.0 

1.0 

2000.0 

50.0 

50.0 

10.0 

40.0 

[((9.5  x  107)/7)exp(- 1150/73]  [12] 
[((4.2  x  107)/73exp(- 1200/73]  [12] 
3.34  x  104exp(- 10300/73  [12] 

0.3  [18] 

4.5  [10] 

1.0  [7] 

15  [21] 

0.1  [18,29] 

0.1  [18,29] 

0.5  [21,23] 

10.78 

2.38 

1320  [7] 

400  [21] 

0.5 

0.5 

0.0171  exp(— 103191/R73 
Afs/Afps 

S"ps  =  exp(— 0.2935?3  +  0.635?2  + 
4.1788?  +  0.3169), 

?  =  (1000/7X10)-  1  [32] 


the  flow  of  reactant  species  through  the  void  spaces.  The  cathode 
overpotential  also  includes  all  forms  of  overpotentials  (activa¬ 
tion,  ohmic  and  concentration).  The  electrolyte  overpotential  is 
due  to  the  resistance  to  the  ion  transport  through  the  electrolyte 
layer.  Although  the  thickness  of  the  cathode  (includes  backing 
and  reaction  zone  layer)  and  electrolyte  are  relatively  thin  in 
anode-supported  SOFC  but  the  potential  loss  in  these  compo¬ 
nents  are  significant  and  hence  cannot  be  neglected. 

The  composition,  morphology  and  thickness  of  the  reaction 
zone  layers  has  a  significant  impact  on  the  cell  performance 
and  a  substantial  amount  of  research  is  being  devoted  in  opti¬ 


mizing  these  layers  [18,21,23,27-29].  The  parameters  which 
are  of  significant  importance  for  enhancing  cell  performance 
are  an  efficient  triple  phase  boundary  (TPB)  area  or  active  area 
for  electrochemical  reaction,  and  high  ionic  and  electronic  con¬ 
ductivities  in  the  ion  and  electron  conducting  particles  [21]. 
According  to  Costamagna  et  al.  [21],  the  largest  active  area  is 
achieved  when  the  dimensions  and  volume  fractions  of  the  ion 
conducting  and  electron  conducting  particles  are  equal.  There¬ 
fore,  in  the  base  case  cell  performance  shown  in  Fig.  3,  the 
dimensions  and  volume  fractions  of  the  ion  conducting  and  elec¬ 
tron  conducting  particles  are  set  to  be  equal. 

Often,  it  is  reported  that  the  contribution  of  concentration 
overpotential  in  anode- supported  SOFC  is  significant  due  to 
thick  anode  [7,10],  but  the  cell  performance  obtained  at  the  base 
case  conditions  reveals  that  the  anode  concentration  overpoten¬ 
tial  is  about  four  orders  of  magnitude  smaller  than  the  anode 
ohmic  overpotential  even  at  higher  current  densities.  This  can 
be  evident  from  Fig.  4.  Moreover,  the  contribution  of  anode 
activation  overpotential  is  found  to  be  negligible.  On  the  other 
hand,  due  to  the  thin  cathode,  the  contribution  of  cathode  con¬ 
centration  overpotential  to  the  total  cell  potential  is  negligible, 
whereas  cathode  activation  overpotential  is  about  two  orders  of 
magnitude  smaller  than  cathode  ohmic  overpotential,  as  shown 
in  Fig.  5.  It  is  clear  from  Figs.  4  and  5  that  the  greatest  con¬ 
tributor  to  the  cell  potential  loss  is  ohmic  overpotential.  More 
specifically,  it  is  the  ionic  conductivity  of  the  ion  conducting 
particles  in  the  reaction  zone  layers,  which  is  contributing  to  the 
overall  ohmic  overpotential,  and  hence  it  needs  to  be  enhanced 
to  improve  the  cell  performance. 

The  effect  of  temperature  on  cell  performance  is  shown  in 
Fig.  6.  The  other  operating  and  design  parameters  are  in  accor¬ 
dance  with  the  base  case  parameters  listed  in  Table  4.  It  can  be 
seen  that  with  the  increase  of  temperature,  the  ideal  or  reversible 
cell  potential  (potential  at  zero  current  density)  decreases, 
whereas  the  actual  cell  potential  increases  with  the  increase  of 
temperature,  in  other  words,  cell  performance  increases  with 
the  increase  of  temperature.  Increasing  the  temperature  not  only 
enhances  the  rate  of  electrochemical  reaction  at  the  reaction  sites 
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Fig.  4.  Anode  overpotentials  at  the  base  case  conditions  given  in  Table  4. 
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but  also  increases  the  electronic  and  ionic  conductivities  of  the 
electron  and  ion  conducting  particles,  respectively,  which  in  turn 
minimizes  the  ohmic  contribution  to  the  total  cell  potential  loss 
and  thereby  enhances  the  cell  performance. 

Fig.  7  shows  the  effect  of  pressure  on  the  performance  of  the 
cell.  The  temperature  and  other  parameters  are  set  as  base  case 
parameters  listed  in  Table  4.  It  can  be  observed  that  increasing  the 
pressure  not  only  increases  the  reversible  cell  potential  but  also 
increases  the  actual  cell  potential.  With  the  increase  of  pressure, 
the  reactant  concentration  at  the  reaction  sites  increases,  which 
in  turn  enhances  the  rate  of  electrochemical  reaction  and  rate 
of  mass  transport  resulting  in  the  minimization  of  anode  and 
cathode  overpotentials  and  hence  better  performance.  However, 
increase  in  pressure  results  in  other  problems  such  as  limitation 
on  material  selection,  gas  sealing  and  mechanical  strength  of  the 
cell  components  [30]. 

The  effect  of  porosity  on  the  performance  of  the  cell  is  shown 
in  Fig.  8.  The  temperature  and  pressure  were  set  at  1073  K  and 
1  atm,  and  all  other  parameters  are  the  same  as  that  of  base  case 


Fig.  7.  Effect  of  pressure  on  cell  performance. 


parameters  given  in  Table  4.  It  can  be  seen  that  increasing  the 
porosity  of  the  cell  components  decreases  the  cell  performance. 
Increasing  the  porosity  increases  the  void  fraction  and  decreases 
the  solid  fraction  of  the  porous  layers  resulting  in  the  reduction 
of  the  active  surface  area  available  for  the  electrochemical  reac¬ 
tion.  Moreover,  the  effective  ionic  and  electron  conductivities  of 
the  porous  layers  decreases  with  the  increase  of  porosity,  which 
results  in  the  increase  of  ohmic  overpotential.  Although  the  con¬ 
centration  overpotential  decreases  with  the  increase  of  porosity 
due  to  the  increased  mass  transport  rates  but  the  cell  performance 
decreases  due  to  increased  ohmic  overpotential  with  porosity. 

Fig.  9  depicts  the  effect  of  tortuosity  on  cell  performance.  The 
temperature  and  pressure  were  again  set  at  1073  K  and  1  atm. 
It  can  be  observed  that  increasing  the  tortuosity  of  the  porous 
layers  decreases  the  performance  of  the  cell.  Increasing  the  tor¬ 
tuosity  of  the  porous  layers  means  increasing  the  tortuous  path, 
which  adds  additional  resistance  to  the  reactant  species  diffusing 
through  the  porous  layers  resulting  in  the  reduction  of  reactant 
concentration  at  the  reaction  sites  and  thereby  decreasing  the  rate 


Fig.  6.  Effect  of  temperature  on  cell  performance. 


Fig.  8.  Effect  of  porosity  on  cell  performance. 
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Fig.  9.  Effect  of  tortuosity  on  cell  performance. 


of  electrochemical  reaction.  In  addition,  the  effective  ionic  and 
electronic  conductivities  decreases  with  the  increase  of  tortuos¬ 
ity  resulting  in  the  increase  of  ohmic  overpotential  and  hence 
cell  performance  decreases  with  the  increase  of  tortuosity. 

The  effect  of  volume  fraction  of  electron  conducting  particles 
of  reaction  zone  layers  on  cell  performance  is  shown  in  Fig.  10. 
Again  the  operating  and  other  design  parameters  are  kept  similar 
to  that  of  base  parameters  listed  in  Table  4.  It  can  be  seen  that 
cell  performance  increases  when  the  volume  fraction  of  electron 
conducting  particles  (0)  in  the  reaction  zone  layers  is  increased 
from  0.4  to  0.5,  and  it  decreases  when  0  is  increased  from  0.5  to 
0.6.  The  increase  and  decrease  in  cell  performance  with  volume 
fraction  of  electron  conducting  particles  demonstrate  the  fact 
that  the  largest  active  area  is  achieved  when  the  dimensions  and 
volume  fractions  of  the  ion  conducting  and  electron  conduct¬ 
ing  particles  are  equal  [21].  When  0  equals  to  0.5,  the  volume 
fraction  of  ion  conducting  particles  (1  —  0)  becomes  equal  to 
the  volume  fraction  of  electron  conducting  particles  resulting  in 


Fig.  10.  Effect  of  volume  fraction  of  electron  conducting  particles  on  cell  per¬ 
formance. 


the  largest  active  area  for  electrochemical  reaction  and  hence 
the  best  performance.  Moreover,  the  effective  ionic  and  elec¬ 
tronic  conductivities  in  the  reaction  zone  layers  are  functions 
of  the  volume  fraction  of  electron  conducting  particles,  increas¬ 
ing  0  beyond  0.5  increases  the  effective  electronic  conductivity 
but  decreases  the  effective  ionic  conductivity,  which  results  in 
increasing  the  ohmic  overpotential  and  thereby  decreasing  the 
cell  performance. 

5.  Conclusions 

A  mathematical  model  of  solid  oxide  fuel  cell  (SOFC)  has 
been  developed,  which  can  predict  the  performance  at  various 
operating  and  design  conditions.  The  model  is  fuel  flexible, 
which  can  predict  the  performance  for  not  only  pure  FF  but 
also  any  reformate  composition  as  fuel.  The  important  feature 
of  this  model  is  the  consideration  of  reaction  zone  layers  as 
finite  volumes,  which  were  often  treated  as  mathematical  sur¬ 
faces  (boundary  conditions)  in  the  existing  models.  In  addition, 
micro  characteristics  of  the  electrodes  are  incorporated  into  the 
model.  The  predicted  performance  is  validated  with  the  exper¬ 
imental  data  found  in  the  literature.  An  excellent  agreement  is 
obtained  between  the  predicted  values  and  the  measured  data 
with  R2 -value  as  high  as  0.988.  Moreover,  the  effect  of  oper¬ 
ating  and  design  conditions  on  the  cell  performance  has  been 
examined.  It  is  found  that  the  most  significant  contribution  to 
the  cell  potential  loss  is  from  the  anode  side  of  the  cell  in  an 
anode-supported  SOFC,  and  anode  concentration  overpotential 
at  base  case  conditions  is  about  four  orders  of  magnitude  smaller 
than  the  anode  ohmic  overpotential  even  at  higher  current  densi¬ 
ties.  Overall,  ohmic  overpotential  is  the  single  largest  contributor 
to  the  total  cell  potential  loss  and  hence  needs  to  be  minimized 
to  enhance  the  cell  performance. 
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